home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / GCF.DEM < prev    next >
Text File  |  1991-04-29  |  733b  |  31 lines

  1. PROGRAM d6r9(input,output,dfile);
  2. (* driver for routine GCF *)
  3. VAR
  4.    a,gammcf,gln,val,x : real;
  5.    i,nval : integer;
  6.    txt : string[25];
  7.    dfile : text;
  8.  
  9. (*$I MODFILE.PAS *)
  10. (*$I GAMMLN.PAS *)
  11.  
  12. (*$I GCF.PAS *)
  13.  
  14. BEGIN
  15.    glopen(dfile,'fncval.dat');
  16.    REPEAT readln(dfile,txt) UNTIL (txt = 'Incomplete Gamma Function');
  17.    readln(dfile,nval);
  18.    writeln(txt);
  19.    writeln('a':4,'x':11,'actual':14,'gcf(a,x)':13,
  20.       'gammln(a)':13,'gln':8);
  21.    FOR i := 1 to nval DO BEGIN
  22.       readln(dfile,a,x,val);
  23.       IF (x >= (a+1.0)) THEN BEGIN
  24.          gcf(a,x,gammcf,gln);
  25.          writeln(a:6:2,x:12:6,(1.0-val):12:6,
  26.             gammcf:12:6,gammln(a):12:6,gln:12:6)
  27.       END
  28.    END;
  29.    close(dfile)
  30. END.
  31.